TIFR/TH/06-28, SINP/TNP/06-25 



o 
o 



Thermodynamics of the PNJL model with nonzero baryon and isospin chemical 

potentials 

Swagato Mukherjee* 
Department of Theoretical Physics, 
Tata Institute, of Fundamental Research, 
Homi Bhabha Road, Mumbai 400 005, India. 

Munshi G. Mustafa^ and Rajarshi Ray-'- 
Theory Division, Saha Institute of Nuclear Physics, 1 /AF, Bidhannagar, Kolkata 700 064, India. 

We have extended the Polyakov-Nambu-Jona-Lasinio (PNJL) model for two degenerate flavours 
to include the isospin chemical potential (/i/). All the diagonal and mixed derivatives of pressure 
with respect to the quark number (proportional to baryon number) chemical potential (/lo) and 
isospin chemical potential upto sixth order have been extracted at = fii = 0. These derivatives 
give the generalized susceptibilities with respect to quark and isospin numbers. Similar estimates 
for the flavour diagonal and off-diagonal susceptibilities are also presented. Comparison to Lattice 
QCD (LQCD) data of some of these susceptibilities for which LQCD data are available, show 
similar temperature dependence, though there are some quantitative deviations above the crossover 
^ ^ temperature. We have also looked at the effects of instanton induced flavour-mixing coming from 

'^'n , the (7a(1) chiral symmetry breaking 't Hooft determinant like term in the NJL part of the model. 

' The diagonal quark number and isospin susceptibilities are completely unaffected. The off-diagonal 

' susceptibilities show significant dependence near the crossover. Finally we present the chemical 

potential dependence of specific heat and speed of sound within the limits of chemical potentials 
• where neither diquarks nor pions can condense. 

> . 

Q>^ PACS numbers: 12.38.Aw, 12.38.Mh, 12.39.-x 
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\^ , I. INTRODUCTION 
O 

' Two most important features of strongly interacting matter at low temperature and chemical potentials are the 

Oh! phenomenon of color charge confinement and chiral symmetry breaking. Ho-wever, -with the increase in temperature 

^ t and/or chemical potential, various phases may appear -with different confining and chiral properties. At present both 

' theoretical and experimental endeavours are under-way to map out the phase diagram of QCD. 

, In the limit of infinite quark mass, the thermal average of the Polyakov-loop can be considered as the order 

^ ' parameter for the confinement-deconfinement transition [1]. Though in presence of dynamical quarks the Polyakov- 

. ^ , loop is not a rigorous order parameter for this transition, it still serves as an indicator of a rapid quark-hadron 

^ ■ crossover. Motivated by this observation, Polyakov-loop based effective theories have been suggested [2-4] to capture 

J-j I the underlying physics of the confinement-deconfinement transition. The essential ingredient of these models is an 

. . . 1 effective potential constructed out of the Polyakov-loop (and its complex conjugate). More recently, the parameters 
in these effective theories have been fixed [5, 6] using the data from Lattice QCD (LQCD) simulations (similar 
comparisons of perturbative effects on Polyakov-loop -with Lattice data above the deconfinement transition -was studied 

With the small quark masses the QCD Lagrangian has a partial global chiral symmetry, which is ho-wever broken 
spontaneously at low temperatures (and hence the absence of chiral partners of low- lying hadrons). This symmetry is 
supposed to be partially restored at higher temperatures and chemical potentials. The chiral condensate is considered 
to be the order parameter in this case. Various effective chiral models exist for the study of physics related to the 
chiral dynamics, e.g. the sigma model [8] and the Nambu-Jona-Lasinio (NJL) model [9, 10]. The parameters of these 
models are fixed from the phenomenology of the hadronic sector. 

Various studies of the QCD inspired models indicate (see e.g. Refs. [11-15]) that at low temperatures there is a 
possibility of first order phase transition for a large baryon chemical potential hb^- This /ib^ is supposed to decrease 



'Electronic address: swagatoOtlf r . res . in 
^Electronic address: munshigolam.mustaf aOsaha. ac . in 
t Electronic address: rajarshi .rayOsaha. ac . in 



2 



with increasing temperature. Thus there is a first order phase transition line starting from (T = 0, ~ M-Bc) the 
axis in the [T,^b) phase diagram which steadily bends towards the (T = T^, ^b = 0) point and may actually 
terminate at a critical end point (CEP) characterized by (T = Te, Mb = i^BE)i which can be detected via enhanced 
critical fluctuations in heavy-ion reactions [16]. The location of this CEP has become a topic of major importance 
in effective model studies (see e.g. Rcf.[17]). For ^b LQCD has a complex determinant which hinders usual 
importance sampling techniques. However recently the CEP was located for the physical [18] and for somewhat larger 
[19] quark masses using the reweighting technique of [20], and for Taylor expansion method in [21]. 

For nonzero isospin chemical potential {^i) models and effective theories [22] find an interesting array of possible 
phases. The most important phenomenon that is supposed to happen is a transition to the pion condensed phase 
close to ~ m^r. This has also been supported by Lattice simulations [23], which does not suffer from the complex 
determinant problem for /x/ ^ and /is = 0. 

In this paper we study some of the thermodynamic properties of strongly interacting matter using the Polyakov loop 
enhanced Nambu-Jona-Lasinio (PNJL) model [25, 26]. In this model one is able to couple the chiral and deconfinement 
order parameters inside a single framework. While the NJL part is supposed to give the correct chiral properties, 
the Polyakov-loop part simulates the deconfinement physics. In fact studies of Polyakov loop coupled to chiral quark 
models have become quite fashionable these days (see e.g. Rcf. [27]). 

The initial motivation to couple Polyakov loop to the NJL model was to understand the coincidence of chiral 
symmetry restoration and deconfinement transitions observed in LQCD simulations [28]. While the NJL part is 
supposed to give the correct chiral properties, the Polyakov-loop part simulates the deconfinement physics. Indeed 
the PNJL model worked well to obtain the "coincidence" of onset of chiral restoration and deconfinement [25, 26]. 
Recently the introduction of the Polyakov loop potential [26, 29] has made it possible to extract estimates of various 
thermodynamic quantities. The pressure, scaled pressure difference at various quark chemical potential /xg (or baryon 
chemical potential hb, where /is = Syito), quark number density and the interaction measure were extracted from 
the PNJL model in Ref. [29] for two quark flavours, and all the quantities compared well with the LQCD data. 
Following this some of us made a comparative study [30] of the quark number susceptibility (QNS) and its higher 
order derivatives with respect to /xg with LQCD data. Here the qualitative features match very well though there are 
some quantitative differences. Very recently the spectral properties of low lying meson states have been studied in 
[31]. 

Encouraged by these results, in this paper we have extended the the PNJL model to incorporate the effects of 
nonzero isospin chemical potential (/x/). The motivation for this is that, it enables one to calculate the isospin number 
susceptibility (INS) and its higher order derivatives with respect to ^q. LQCD data on these quantities are also 
available [32] . Thus comparing the results of PNJL for these quantities with that for the LQCD data will provide an 
opportunity to perform some stringent tests on the PNJL model. 

Moreover, once both the QNS and the INS are known one can proceed further to compute the flavour diagonal and 
off-diagonal susceptibilities separately. Since the 2-nd order flavour off-diagonal susceptibility measures the correlation 
among "up" (u) and "down" (d) flavours [33] , this quantity provide a direct understanding to the extent in which the 
PNJL model captures the underlying physics of QCD. 

In our attempt to have a closer look at the u-d flavour correlation within the PNJL model, we have modified the 
NJL part of the PNJL model by using the NJL Lagrangian proposed in [34] . This Lagrangian has a term that can be 
interpreted as an interaction induced by instantons and reflects the [/yi(l)- anomaly of QCD. It has the structure of a 
't-Hooft determinant in the flavour space, leading to flavour-mixing. By adjusting the relative strength of this term 
one can explicitly control the amount of flavour-mixing in the NJL sector. This modified NJL Lagrangian reduces 
to the standard NJL Lagrangian [9, 10] in some particular limit. This modification of the PNJL model have allowed 
us to study the effects of such flavour-mixing on various susceptibilities, specially on the 2-nd order off-diagonal one 
which measures the u-d flavour correlation. 

Investigation of the flavour-mixing effects brings us to an important issue regarding the NJL-type models. Within 
the framework of an NJL model it has been found [35] that for jij = 0, in the T — fiQ plane, there is a single first 
order phase transition line (which ends at a critical endpoint) at low temperatures. But for /i/ 7^ this single line 
separates into two first order phase transition lines because of the different behaviour of the u and d quark condensates 
[36]. Thus there is a possibility of having two critical end-points in the QCD phase diagram [36]. This has also been 
observed in Random Matrix models [37] , in ladder QCD models [38] as well as in hadron resonance gas models [39] . 
It was then argued in Ref. [34, 40] that the flavour-mixing through the instanton effects [41, 42] may wipe out this 
splitting. Later studies found that the splitting is considerable when /i/ is large [39, 43] or fiB is large [44]. We shall 
restrict ourselves only to small chemical potentials and calculate the susceptibilities with the modified PNJL model 
for different amount of flavour-mixing. Comparing these with LQCD data may give us some idea about actual amount 
of the flavour-mixing that is favoured by the LQCD simulations. 

Our next objective is to study the specific heat at constant volume (Cy) and speed of sound («<,) of strongly 
interacting systems. These two quantities are of major importance for heavy-ion collision experiments. While Cy 
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is related to the event-by-event temperature fluctuations [45] and mean transverse momentum fluctuations [46] in 
heavy- ion collisions, the quantity Vg controls the expansion rate of the fireball produced in such collisions and hence 
an important input parameter for the hydrodynamic studies [47-50] . The temperature dependence of these quantities 
were reported earlier in Ref. [30]. For the sake of completeness, in this paper we have also studied the quark number 
and isospin chemical potential dependence of Cy and v^. 

The plan of this paper is as follows. In Section II, we will present our formalisms. First, we will briefly discuss 
the extended PNJL model which wc are going to use. Next, in the same section, formalisms regarding the Taylor 
expansion of pressure (with respect to and fii) and formulae for specific heat Cy and speed of sound Vg will be 
given. In Section III we will present our results and compare some of those with the available LQCD data. Finally, 
we conclude with a discussion in Section IV. Detail mathematical expressions regarding the model can be found in 
Appendix A. 



II. FORMALISM 



A. PNJL Model 



The PNJL model at nonzero temperature T and quark number chemical potential /io was introduced in Ref. [26, 29]. 
Here we extend it to include the isospin chemical potential We have introduced separate chemical potentials fiu 
and Hd for the "up" and "down" quark flavours respectively in the NJL model following Ref. [34, 35]. To further 
extend it to include the Polyakov loop dynamics we have followed the parameterization of the PNJL model used in 
Ref. [29]. Wc start with the final form of the mean field thermodynamic potential per unit volume that we have 
obtained. It is given by (further details about the model can be found in Appendix A), 

2tJ -p^|ln[l + 3(a> + $e-(^/-«)/^)e-(^/-w)/^+e^3(B,-M/)/T 



f=u,d 



+ ln[l + 3($ + $e-(^^+^^)/^)e-(^^+«)/^ + e-3(S/+w)/T]|_ ^ 6 /'^i?/^(A2-p2) . (i) 



Here for the two flavours the respective quark condensates are given by (t„ uu > and ad =< dd > ^ and the 
respective chemical potentials are fiu and fid- Note that fio = {fiu + fJ-d)/2 and fii ~ {fiu — Md)/2. The quasi-particle 



energies are Eu,d = yP^ + "^u d' where mu,d = rno — 4:GiUu,d — '^G2<Jd,u are the constituent quark masses and mo 

is the current quark mass (we assume flavour degeneracy). Gi and G2 are the effective coupling strengths of a local, 
chiral symmetric four-point interaction. A is the 3-momentum cutoff in the NJL model. U ($, i>, T) is the effective 
potential for the mean values of the traced Polyakov-loop $ and its conjugate and T is the temperature. The 
functional form of the potential is. 



U ($,«•, T) 
2^4 



2 6 
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(2) 



with 



&2(T) = ao + ai (I) +a3(|) ■ (3) 

The coefficients and hi were fitted from LQCD data of pure gauge theory. The parameter Tq is precisely the 
transition temperature for this theory, and as indicated by LQCD data its value was chosen to be 270 MeV [51-53]. 
With the coupling to NJL model the transition doesn't remain first order. In this case from the peak in d^/dT the 
transition (or crossover) temperature comes around 227 McV. 

Before we move further we note some important features of this model: 



^ Here we deviate from the convention of defining the sigma condensates from those of Refs. [29, 30]. 
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• Since the gluons in this model are contained only in a static background field, the model would be suitable to 
study the physics below T = 2.5Tc. Above this temperature the transverse degrees of freedom become important 
[54]. 

• In general, pion condensation takes place in NJL models for /i/ > Also there is a chiral transition 
for fiQ ~ 340 MeV above which diquark physics become important. For simplicity we neglect both the pion 
condensation and diquarks ^ and so restrict our analysis to < 70 MeV and /ig < 200 MeV. 

• As discussed in the appendix, for G2 = the full symmetry of the Lagrangian in the chiral limit (mo = 0) is 
SUv{2) X SUa{'2) X C/y(l) x [/^(l). The coefficient G2 is interpreted as inducing instanton effects as it breaks 
the 1/^(1) symmetry explicitly by mixing the quark flavours. By using a parameterization Gi = (1 — a)Go and 
G2 = (xGq (following Ref. [34]), one can tune the amount of instanton induced flavour mixing by varying a. For 
a = there is no instanton induced flavour mixing, and for a = 1 the mixing becomes maximal. We shall look 
into the effects of this mixing in the susceptibilities. 

The form of the NJL part in Eqn. (1) is a generalization of the standard NJL model, which we get when Gi = G2, 
and jiu = fJ'd- In fact the potential in Eqn. (1) becomes exactly the same as that of Ref. [29, 30] if we use a = 0.5 
and put Go equal to half the four-point coupling G in those references. We shall use this value for Go in this 
work. 

• For the NJL sector without coupling to the Polyakov loop (i.e. setting $ = $ = 1) one can easily see that 
the expression for f2 in Eqn. (1), is invariant under the transformations fiu —fJ-u "and/or" fid — > —fid- This 
implies that the physics along the directions of /xo = and = at any given temperature are equivalent. 
However inclusion of the Polyakov loop turns off this symmetry. Now f2 is invariant only under the simultaneous 
transformation —fJ-u "and" /i^ —fJ-d- This is a manifestation of the CP symmetry which implies that 
n is symmetric only under the simultaneous transformation $ — > $ and fi^^d ^ ^Mu.d and vice-verse. Thus 
coefficients of $ and $ are found to be equal when fiQ = 0, and different when fii ~ 0. In the T — fij plane 
we shall have $ = $, and everywhere else $ ^ $. This is reminiscent of the complex fermion determinant for 
nonzero /xq. This will be seen to have important consequences for the extraction of susceptibilities. 

• On the other hand the quark condensates cJu and ad are equal to each other whenever /io = or /i/ = [see, e.g., 
Eqn. (A18)]in the NJL as well as PNJL model. This can be seen by inspecting the thermodynamic potential 
and remembering that we are using Gi + G2 = Gq, and also the fact that for /io = 0, $ = $. Now Gi and G2 
are only coupled to the Uu and ad- It is clear from Eqn. (1) that whenever cr„ = ad, the couplings Gi and G2 
come in the combination Gi -I- G2 = Go = constant. This means that the physics is completely independent of 
these couplings whenever either /io = or /*/ = 0. 

B. Taylor expansion of Pressure 

The pressure as a function of temperature T, quark chemical potential /ig and isospin chemical potential /// is given 



P(r,/io,/i/) = -f}(T,/io,/i/) . (4) 

Following usual thermodynamical relations one can show that the first derivative of pressure with respect to /to gives 
the quark number density. The second derivative is the quark number susceptibility. In LQCD since usual Monte 
Carlo importance sampling fails for nonzero /ioi the QNS and higher order derivatives computed at /io = can be 
used as Taylor expansion coefficients to extract chemical potential dependence of pressure. 

Given the thermodynamic potential f2, our job is to minimize it with respect to the fields ct„, ad, and $, using 
the following set of equations, 

dn an an 

9^^°' 9$^°' m^'- 



^ Very recently diquarks have been discussed in Ref. [55] and pion condensation in Ref. [56]. 
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The values of the fields so obtained can then be used to evaluate all the thermodynamic quantities in mean- 
field approximation. The cross-over temperature for /xq = /x/ = was obtained in Ref. [30] and was found to be 
Tc = 227 MeV. The field values obtained from Eqn. (5) arc then put back into fl to obtain pressure from (4). We can 
then expand the scaled pressure at a given temperature in a Taylor series for the two chemical potentials /io and /i/, 



rp4 



n=0 j=0 •' ^ ■' ' 



where, 



1 d-^{P{T,p.^,^Li)lT^ 



(6) 



(7) 



The n = odd terms vanish due to CP symmetry. Even for the n = even terms, due to flavour degeneracy all the 
coefficients with j and k both odd vanish identically In this work we evaluate all the 10 nonzero coefficients 
(including the pressure at /io = /ij = 0) upto order ti = 6. Some of these coefficients have already been measured on 
the LQCD [24, 32]. In our earlier work [30], we compared the 4 coefficients for /i/ = with those of the LQCD data 
using improved actions [32]. Here we shall be able to compare 3 more coefficients with LQCD data and also predict 
the behaviour of the other 3 coefficients. 

Let us now identify the coefficients which we shall compare with the LQCD data. The first set is given by. 



c„(T) 



\_ 9" {P{T,p.^)lT^ 



Mo=0 



These coefficients were already computed upto 8-th order and compared to LQCD data to 6-th order in 
new set of coefficients to be compared with the LQCD data upto n = 6 are, 



1 a" (p(r,/io,/i/)/r4 



.("-2)2 



n> I. 



(8) 



The 



(9) 



fiQ=a,fj.i=o 



The remaining coefficients we obtain are 04^, and Cg^. 

To complete the comparison with the LQCD data we have looked at the flavour diagonal (c"") and flavour off- 
diagonal {c!^'^) susceptibilities defined as. 



„nO , „("-2)2 
L„ -r On 



and 



" 4 ' " 4 

The 2-nd order flavour diagonal and off-diagonal susceptibilities are given by, 

XuuiT, = 0, = 0) d'^P{T, Md) 



(10) 



T2 

XudiT, flu = 0,/id 



5m« 

0) d^PiT,t,^,fXd) 



and 



j<2 



P.u=tJ.d=0 



2c"''. 



In this work we have computed all the coefficients using the following method. First the pressure is obtained as 
a function of /io and fij for each value of T, and then fitted to a sixth order polynomial in /io and /i/. The quark 
number susceptibility, isospin number susceptibility and all other higher order derivatives are then obtained from the 
coefficients of the polynomial extracted from the fit. In the fits we have used only the even order terms. 

C. Specific lieat and speed of sound 

Given the thermodynamic potential fi, the energy density e is obtained from the relation. 



d{n/T) 



dT 



= -T — 
dT 



+ n 



(11) 
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The rate of change of energy density e with temperature at constant volume is the specific heat Cv which is given as, 



T 



9T2 



(12) 



For a continuous phase transition one expects a divergence in Cy, which, as discussed ear her, wiU translate into 
highly enhanced transverse momentum fluctuations or highly suppressed temperature fluctuations if the dynamics in 
relativistic heavy-ion collisions is such that the system passes close to the critical end point (CEP) in the T — /ig 
plane. 

The square of velocity of sound at constant entropy S is given by. 



dP 

a7 



dP 

of 



df 



an 

df 



T 



gf2 



(13) 



Since the denominator is nothing but the Cy, a divergence in specific heat would mean the velocity of sound going 
to zero at the CEP. 

Given the relations Eqn. (12) and Eqn. (13), we first obtain the f2(T, /xq = 0) from the PNJL model. We then 
obtain the derivatives using the standard finite difference method. To get points close enough we have used cubic 
spline interpolations. This procedure has been repeated for various values of and 

III. RESULTS 

A. Taylor expansion of Pressure 

As discussed in section II B, we extract the Taylor expansion coefficients by fitting the pressure as a function of fiQ 
and Hi at each temperature. Data for pressure was obtained in the range < /io < 50 MeV and < /z/ < 50 MeV at 
all the temperatures. Spacing between consecutive data was kept at 0.1 MeV. We obtain all possible coefficients upto 
6*'' order using gnuplot program. The least-squares of all the fits came out to be 10"^'' or less. This method was 
already used in our earlier work [30] where we checked the reliability of such fits. Here again we have reproduced all 
those coefficients satisfactorily. We shall first discuss the results with the standard flavour mixing in the NJL model 
parameterization (i.e. with Gi = 6*2), and then discuss the results for minimal (G2 = 0) and maximal (Gi = 0) flavour 
mixing. 



Gi ~ G2 



We start by presenting our results for the PNJL model with the standard NJL Lagrangian, i.e., Gi = G2 = Go/2. 
Note that this is the case studied in the PNJL models of Refs. [26, 29, 30], but without the isospin chemical potential. 

We present the QNS, INS and their higher order derivatives with respect to /io in Fig. 1. We have also plotted 
the LQCD data from Ref. [32] for quantitative comparison. At the second order of Taylor expansion i.e. n = 2 we 
find (also observed earlier in [30]) that the QNS C2 compares well with the LQCD data. On the other hand, the INS 
C2 quickly reaches its ideal gas value above Tc (around 2Tc) in our model calculations, whereas the LQCD value are 
lower and matches with the value of C2. Note that in the present form of the model the Polyakov loop itself rises a 
little above 1 and saturates. This leads to the INS to rise slightly above 1 at high temperatures. At the 4*'* order we 
see that the values of C4 (also observed in [30]) in the PNJL model matches closely with those of LQCD data for upto 
T ~ 1.05Tc and deviates significantly thereafter. The coefficient is close to the LQCD data for the full range of T 
upto 2Tc. 

Earlier expectation [30, 57] was that, the mean field analysis may not be sufficient and hence the higher order 
coefficient C4 in the PNJL model shows significant departure from lattice results. This should have also meant that 
the INS C2 should be more closer to LQCD data than C4. However, our results show that the INS C2 is significantly 
different from the LQCD data above Tc, but c{ is quite consistent. Further, we sec from Fig. 1 that both the 6*'' 



^ see http://www.gnuplot.info/ 
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FIG. 2: Left panel: $ (solid lines) decreases and $ (dotted lines) increases as a function of fio/T {fii — 0). Right panel: $ 
(solid lines) and $ (dotted lines) are equal and almost constant as a function of fii/T {no = 0). 



order coefficients cg and Cg are quite consistent with the LQCD results. We now give a qualitative explanation for 
the PNJL results and try to understand the behaviour of the coefficients above T^. 

We pointed out in Section II A that in the thermodynamic potential Eqn. (1), the Polyakov loop couples to /io ^^nd 
its conjugate couples to — yUo due to CP symmetry. As observed in SU{N) matrix model [58] and also in the PNJL 
model [29, 30], this difference in coupling leads to splitting of the Polyakov loop and its conjugate for any nonzero fio- 
Thus even at high temperatures when the Polyakov loop is close to 1, it decreases with increasing /io and its conjugate 
increases (see left panel of Fig. 2). This means that the /io dependence of pressure is not the same as that for an 
ideal gas. Hence the coefficients C2 and C4 arc both quite off from their respective ideal gas values. Also note that 
though ce is close enough, it is still distinctly different from zero. On the other hand for /io = 0, the Polyakov loop 
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as well as its conjugate couples to both the ^/ and — yUj. They are, thus, equal (see right panel of Fig. 2), and also 
found to be almost constant for small So the temperature dependence of the INS and its ^/ derivatives should 
reach the ideal gas behaviour above Tc. For the coefficients which are mixed derivatives of and the behaviour 
should be somewhere in between. And indeed we see that C2, c\ and Cg in Fig. 1 are quite close to their respective 
ideal gas values above Tc. Thus the LQCD results that show almost equal values of QNS and INS, indicate that the 
splitting between the Polyakov Loop and its conjugate in the /i/ = direction for T > l.5Tc is almost negligible (also 
supported by pQCD). This splitting was taken to be absolutely zero in the recent report with the PNJL model in 
Ref. [55]. 




0.6 0.8 1 1.2 1.4 1.6 1.8 2 0.6 0.8 1 1.2 1.4 1.6 1.8 2 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

T-To m„ T/T, 



FIG. 3: The flavour diagonal {upper row) and flavour off-diagonal {lower row) susceptibilities for n = 2, 4 and 6 as functions of 
T /Tc. Symbols are LQCD data [32]. The arrows on the right indicate the respective ideal gas values. 

In order to investigate these discrepancies between the results form the PNJL model and the LQCD data more 
closely, we have also calculated the flavour diagonal {c^) and off-diagonal {djf) susceptibilities, defined in Section 
II B, upto 6-th order. These are shown in Fig. 3. Except for Cj", all the other LQCD results for flavour diagonal 
susceptibilities are close to their respective ideal gas values from around 1.2Tc onwards. The PNJL model values for the 
diagonal coefficient C2 " seem to be more or less consistent with the LQCD data. The most striking discrepancy with 
the LQCD data shows up in the 2-nd order flavour off-diagonal susceptibility Cj''. As discussed earlier, c^'^ signifies 
the mixing of u and d quarks through the contribution of the two disconnected u and d quark loops. While the LQCD 
data shows that this kind of correlation between the u-d flavours are almost zero just away form Tc, the PNJL model 
results remains significant even upto 2Tc. The negativity of 02'^ (see Fig. 3) indicates that the dominant correlation 
is between u quarks and d anti-quarks and vice-verse, i.e., pion-like. Hence putting in the dynamical pion condensate 
may throw some light on this issue. Also addition of any new couplings (e.g. as shown for the isoscalar-vector and 
isovector- vector couplings for NJL model in Ref. [59]) may have important consequences for these suscpetibilities. 

Again from Fig. 3, the 4-th order diagonal (c™) as well as the off-diagonal (04'') coefficients show a behaviour 
similar to c|°. Whereas the LQCD data reaches the ideal gas value above Tc, the PNJL values are quite distinctly 
separated. Finally, at the Q*^ order the behaviour for both diagonal and off-diagonal coefficients in the PNJL model 
and LQCD are quite consistent. 

We now present the temperature dependence of the remaining nonzero coefficients (Fig. 4). is the 4-th order 
diagonal coefficient in the isospin direction. In contrast to we see that approaches the ideal gas value quite 
fast above Tc. The behaviour of Cg^ is quite similar to its counterpart c^'^ . This is in accordance to the expectation, 
as discussed earlier. Same is true for the coefficient c|"'. 
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FIG. 4: cT and cT and Cg as functions of T/Tc. Arrows on the right indicate the respective ideal gas values. 




T/Tc T/Tc T/Tc 

FIG. 5: Left panel: £2^* is independent of a. Middle and Right panel: Dependence of some off-diagonal coefficients on the 
flavour mixing parameter a. Symbols are Lattice data [32]. Arrows on the right indicate the corresponding ideal gas values. 

2. Gi / Ga 

Since the PNJL model has problem in reproducing the LQCD data for c!^'^ which is a measure of the flavour-flavour 
correlation, it is interesting to have a closer look at the effect of flavour-mixing on different susceptibilities. As 
discussed earlier, the parameterization Gi = (1 — a)Go and G2 = (xGo enables one to tune the instanton induced 
flavour mixing by varying the value of a between and 1. Here we discuss the two extreme cases of a = 1 (maximal 
mixing) and a = (zero mixing). We have re-calculated all c„ and c^, upto n = 6, for a = 0, and, 1. We found that 
all the diagonal coefficients, including and C4 whose behaviour are the most drastically different in the PNJL model 
and in LQCD, are independent of the values of a. As a consequence, the 2-nd order flavour off-diagonal susceptibility 
[cj'^ = — C2)/4] is also unaffected by the instanton induced flavour- mixing effects. 

The above fact can be understood from the following reasoning. We mentioned in section II A that the quark 
condensates (Ju and cr^j are equal to each other for either of the cases — and /i/ = 0. This is clear from Eqn. A18. 
Now Gi and G2 couple only to the <7u and a^,- So for (t„ = ad, we only get the combination Gi -|- G2 = Go, which is 
a constant. Thus none of the physics in the /zq = and ^/ = directions depend on the value of a, implying that 
the diagonal derivatives in these two directions will also be independent of a. 

However, the mixed derivatives can have dependence on a. This is because the values of cr„ and ad can be different 
when both fiQ and /i/ are together nonzero. This was seen in Ref. [34] for the normal NJL model. But those authors 
also found that there is a critical value of ac « 0.11 above which the condensates (7„ and ad become equal even for 
both HQ and being nonzero. Here, for the PNJL model we have found that all the mixed derivatives upto 6*'* order 
are exactly equal for the two cases a = 0.5 (standard mixing used in NJL and PNJL models) and a ~ 1 (maximal 
mixing) which is in accordance to the results of the above reference. We hope to obtain the value for ac for the PNJL 
model in future. For a = all the off-diagonal coefficients were found to differ from those at a = 0.5. 

The left-most panel in Fig. 5 shows the independence of €2"^ on a. The rest figures show one representative 
coefficient each for n = 4, and 6. As can be seen, the instanton effects quite significantly suppress the temperature 
variation of these coefficients near Tc. Also it can be observed from Fig. 5, that the LQCD data favours larger amount 
of instanton induced flavour-mixing. 
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B. Dependence of Cv and v1 on /io and /i/ 
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FIG. 6: Cv as a function of T/Tc- Left panel shows the variation with ^o; Right panel shows the variation with /i/. Arrows on 
the right indicate the ideal gas value for jj.o — fJ-i — 0. 
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FIG. 7: as a function of T/Tc- Left panel shows the variation with ^o; Right panel shows the variation with fij. Arrows on 
the right indicate the ideal gas value for no — fii — 0. 

Here we present the chemical potential dependence of specific heat Cy and the speed of sound Vg . The range of the 
three representative values of fio and /i/ are such that neither the diquark physics nor the pion condensation becomes 
important. In the ideal gas limit the expression for Cv is as a function of temperature T and either of the chemical 
potentials fiQ or fij is given by, Cv/T^ = (747r^/15) + 6{fil j/T^). Thus, for large temperatures and not so large 
chemical potentials, it can be expected that the Cy is more or less independent of Hoj^s. This is borne out in the 
PNJL model as seen in Fig. 6. At low temperatures however, there can be non-trivial contribution from chemical 
potential. As illustrated in Fig. 6, the low temperature behaviour is away from ideal gas, but there is significant 
difference in the values of Cy as a function of fif). In the range of fij considered, even for T < Tc there seems to 
no significant isospin effects. Another interesting feature is that as a function of /ip, the peak of Cy which appears 
at Tc shifts towards lower temperatures. This signifies that the transition temperature may decrease and also the 
nature of transition may change as the chemical potentials increase. A decrease of Tc with increasing /xg and is 
consistent with what have been found on the Lattice [60, 61]. We hope to address this issue through the analysis of 
chiral susceptibility in a future publication. 

The speed of sound in the ideal gas limit is the same -y/S for any given temperature and chemical potential. As 
shown in Fig. 7 the v1 for different //q and merges towards the ideal gas value at large temperatures. However, 
even above T^ there is significant increase in for increase in /ig. So for nonzero quark matter density the speed 
of sound is higher near Tc and this may have important contribution to thcrmalization of the matter created in 
relativistic heavy-ion collision experiments. Again there seems to be negligible isospin dependence of in the range 



11 



of temperatures studied. From Fig. 7 we note that in the PNJL model even with /ip as large as O.STc, the never 
reaches a value as large as 0.2 near or below T = Tc which was used in [62] to describe the rapidity spectra. 

IV. DISCUSSIONS AND SUMMARY 

We have extended the PNJL model of Rcf. [29] by the introduction of isospin chemical potential. Using this we 
have studied the behaviour of strongly interacting matter with two degenerate quark flavours in the phase space of 
T, /io and /i/, for small values of the chemical potentials. We have extracted 10 coefficients of Taylor expansion of 
pressure in the two chemical potentials upto 6-th order. Some of these coefficients were compared with available 
LQCD data. The quark number susceptibility and isospin susceptibility show order parameter-like behaviour. A 
quantitative comparison shows that the quark number susceptibility reaches about 85% of its ideal gas value upto 
temperature of about 2Tc, consistent with LQCD results. However, the isospin susceptibility reaches its ideal gas 
value by this temperature. This is in contrast to LQCD results where both the susceptibilities are almost equal from 
around 1.2Tc onward. Similarly, the higher order derivatives for approach the the ideal gas behaviour much faster 
compared to those for /^o- In contrast, though both the QNS and INS in LQCD deviate from their ideal gas values, 
the higher order derivatives reach their ideal gas limit quickly. The values of the mixed derivatives in the PNJL model 
shows a behaviour somewhat in between. On the Lattice however, the mixed derivatives are almost zero (i.e., the 
ideal gas value) above Tc. 

Thus some of the coefficients in the PNJL model differ from the LQCD data and one could hold the mean field 
analysis responsible for this departure. But if this argument were true then the higher order derivatives obtained in 
the PNJL model should depart from the LQCD data more than the lower order coefficients, which is not the case. As 
against this expectation, we have found a very nice pattern in the PNJL results which can be understood in terms 
of the behaviour of the Polyakov loop. The dependence of the Polyakov loop and its conjugate on temperature and 
the chemical potentials is extremely important. For fij = they have different values when fiQ is varied. This makes 
all the coefficients which are derivatives of pressure with respect to fiQ alone, to deviate from the ideal gas behaviour. 
For jiQ = however the Polyakov loop and its conjugate are equal and hence both reach the ideal gas value above Tc- 
Thus the coefficients which are derivatives of pressure with respect to fij alone, all reach their respective ideal gas 
values above Tc. The mixed derivatives are found to be somewhere in between. Nonetheless, we hope to look into the 
effects of fluctuations in future. 

In order to have a closer look at the discrepancy between the PNJL results and LQCD data, we have also calculated 
the flavour diagonal and flavour off-diagonal susceptibilities upto 6-th order. We have found that the 2-nd order flavour 
off-diagonal susceptibility, which indicates the correlation among "up" and "down" quarks, is significantly away from 
zero even upto T = 2Tc. On the other hand, LQCD results [63] show that correlation among the flavours in the 2-nd 
order off-diagonal susceptibility is largely governed by the interaction of the quarks with the gauge flelds and is almost 
independent of the presence of the quarks loops. This motivated us to study the instanton induced flavour-mixing 
effects within the framework of the PNJL model. Unfortunately, we found no effect of flavour-mixing on any diagonal 
QNS and INS, and hence on the 2-nd order flavour off-diagonal susceptibility. We speculate one possibility to reconcile 
PNJL and LQCD data, that is to keep the pion condensate as a dynamical variable and perform the calculations. In 
fact there are indications [64] that at zero temperature and in the chiral limit, pion condensation can be catalysed by 
an external chromomagnetic field. We hope to present these results in future. 

On the other hand, flavour-mixing effects on the mixed susceptibilities of quark and isospin chemical potentials 
indicate that large flavour-mixing is favoured by the LQCD data. This may have important consequences [34] on the 
phase diagram of the NJL model at low temperature and large baryon chemical potential. 

Apart from the possible improvements for the Polyakov loop potential, inclusion of pionic and diquark fluctuations, 
etc. we also intend to include terms in the NJL part with six point couplings to take proper account of the quark 
number fluctuations in the low temperature phase. 

We have also investigated chemical potential dependence of spcciflc heat and speed of sound. The speciflc heat 
sort of becomes independent of chemical potential just above Tc. Below Tc there is some significant effect from both 
the chemical potentials /io and /i/. Consistent with LQCD findings [60, 61], the peak in specific heat towards lower 
temperatures with increasing chemical potentials indicating a decrease in the transition temperature. We plan to 
make a more detailed investigation of the location of the phase boundary. The speed of sound, on the other hand, 
increases with the increase of either of the chemical potentials in almost the whole range of temperatures. But this 
dependence become milder as one goes to higher temperatures. Thus with a proper implementation of the PNJL 
equation of state into the hydrodynamic studies of elliptic fiow, one may able to make some estimates of both the 
temperature and densities reached in the heavy- ion collision experiments. 
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APPENDIX A 

Here we present the complete details of the PNJL model used in our work. First we discuss the NJL model in the 
complete space of temperature T and the "up" and "down" flavour chemical potentials ^„ and (or equivalently 
the quark chemical potential /io and isospin chemical potential /i/) [34, 35]. Then we extend it to couple with the 
Polyakov Loop. 

1. The NJL model 



The NJL model Lagrangian for two flavours can be written as [34, 35] : 



C 

Co 



Co = G 



C0 + C1+ C2 
2 



where, 



[Gi] = [G2] - [energy]- 



m = diag(mu, md). 



(Ala) 
(Alb) 

(Ale) 
(Aid) 



(A2) 



We shall assume flavour degeneracy to„ = rud ~ m^. For toq = the symmetries of the different parts of the 
Lagrangian Al are: 



Co 
Ci 
C2 



SUvi2) X SUAi2) X Uv{l) X [/^(l) 
SUv{2) X SUa{2) X Uv{l) x C/^(l) 
SUv{2) X SUa{2) xUv{l) 



(A3a) 
(A3b) 
(A3c) 



£2 has the structure of a 't-Hooft determinant, det [g(l + j^Jql+dct [q{l — 75)9] [10], and breaks J7a(1) axial symmetry. 
This interaction can be interpreted as induced by instantons and reflects the [/^(l)-anomaly of QCD. 

We are interested in the properties of this Lagrangian at nonzero temperatures T and chemical potentials /i„ and 
fid- Equivalently, one can also use the quark number chemical potential ^0 = (Mm + t^d)/'^ and the isospin chemical 
potential fij ~ (/i„ — fid)/'2. In the mean field approximation we consider the two quark condensates cr„ = (uu) and 
ad = {ddj. The pion condensate tt is assumed to be zero (which is true in the NJL model for fij < m^jl^ Then the 
thermodynamic potential is obtained as. 



r2(T,^„,/id) = ^ r2o(T,^/;m/) + 2Gi (cr^ +CT^) +4G2cr„crd , 



-2iV< 

-2NrT 



J [Ztt ) 



(Stt) 
(2^ 



In 



1 + -A'/)/T 



(A4a) 



(A4b) 



where, the energy Ef and constituent quark mass m/ is given by. 



irif — mo — 4Gicr/ — 4G2(7/' , 



f^f^{u,d} 



(A5) 
(A6) 
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Finding the stationary points of the thermodynamic potential with respect to (t„ and ad, i.e., solving the coupled 
equations dfl/dau = and dVL/dud = 0, one gets the gap equations, 

af = -2NJ^'^[e{k^-p^)-n{E^)-n{E^)\, f = u,d (A7a) 
J (2tt) J^f 

n(Ef) = ^ -, and niEf) = ^ -. (A7b) 

^ ^' l + e^p{Ef - ^if)' ^ ^' l + exp(£;/+/i/) ^ ' 

The constituent mass nif for one flavour depends in general on both the condensates [see Eqn. (A6)] and therefore 
the two flavours are coupled. Chiral symmetry {SUa{'2)) is broken spontaneously for cr/ ^ 0. 
Let us now make the parameterization 

Gi = {l-a)G„, and Ga = aGo (A8) 

with a fixed value of Gq. Tuning the value of a one can control the flavour mixing in the Lagrangian. We consider 
some of the cases below. 

1. g = : This implies G2 = i.e.the [^4(1) symmetry breaking term £2 drops out and hence C has no flavour 
mixing. 

2. a = 1 : Here Gi ~ 0, and thus £2 completely dominates the coupling. The flavour mixing in C is thus "maximal" . 



3. a — 1/2 : In this case we have Gi = G2 = Go/2. So the Lagrangian £ = £0 + Go {iJH')'^ + (V'*75'^V')^ 



is the 



standard NJL model [9]. Here also the C^a(1) symmetry is broken which is commensurate with the fact that in 
nature the i] particle is much heavier than the tt's 

2. Extension to PNJL 

Our aim is to extend the PNJL model introduced in Ref. .[26, 29] to include isospin chemical potential. To achieve 

this we now include the Polyakov loop and its effective potential to the NJL model described above. The Lagrangian 
becomes, 

CPNJL = Co + Ci+C2-U{^AmAl T) ( A9) 
The only part of the NJL sector that is modified is £0 which now becomes, 

Co = '^l'{ip-m)^lj , (AlO) 

where 

^?^ = 9^-^v4^ A'^^Sf^oA", A''{x) = gAii{x)Xa/2. (All) 

A'^{x) are SU{3) gauge fields and Xa are Gell-Mann matrices. 

U {^[A], ^\A],T) is the effective potential expressed in terms of the traced (over color) Polyakov loop (with periodic 
boundary conditions) and its charge conjugate — 



'0 

i I dr ^4(2;, r) 







/?=i Ai^zAo. (A12) 



We shall be working in the mean field limit. For simplicity of notation we shall use <I> and <& as their respective mean 
fields. $ is the order parameter for deconfinement transition. In the absence of quarks <& = 4> and deconfincment is 
associated with the spontaneous breaking of the Z{3) symmetry. Conforming to this symmetry and parameterizing 
the LQCD Monte Carlo data one can write down an effective potential for $ and Following Ref. [29], we write 

. _fe2),4 - I + 4=) + I . (A13.) 



b,{T) ^ ao + a,[f +a,if + ^3 if ] . (A13b) 
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At low temperature lA has a single minimum at $ = 0, while at high temperatures it develops a second one which turns 
into the absolute minimum above a critical temperature Tg. $ and $ will be treated as independent classical fields. 
The mean field analysis of the NJL part of the model proceeds in exactly the same way as in the previous case. Using 
CT„ and CTd as the independent quark condensates (and neglecting tt) one gets the expression for the thermodynamic 
potential, 



f=u,d 



(Al4a) 



^2Nr I -^Ef0{A^ ^p^)~2T I ^^Tr^ln 



{2nY 
_2r / ^^Tr.ln 



l + Le-(^/-^^)/^ 



1 + Lte"(-^^+^^^/'^ 



(A14b) 



where Eu,d = Y"i^,d + and m„,d = mo - 4:Gi<7u,d ~ 4G2crd,„. 

Note that with a = 0.5 and /i„ = /id = /i and if for the coupling G and condensate a of Ref. [29] one uses G = 2Go 
and (T = au + (Td, then the thermodynamic potentials here and in Ref. [29] are exactly equal. 
Using the identity TrlnA" = IndetX one can write for a given flavour /, 



l + ite-(B/+M/)/T 



Indet 1 + Le-'^^'->''>'-' +lndet 

= In [l + 3 ($ + ^e-^^z-A'/)/^) e-^E,-^,)/T ^ g-3(iJ,-„)/T 
+ In 1 + 3 + $e-(^/+w)/T^^ f,-(Ef+f,f)/T j_ g-3(£;/+M/)/T" 

This gives us the final form of the thermodynamic potential as, 

f=u,d 

J (27rj 



-2T 
-2T 



(2^ 

{27Tf 



In 



In 



(A15) 



(Al6a) 



(A16b) 



From this thermodynamic potential the equations of motion for the mean fields cr„, drf, <f> and $ are derived through, 
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This coupled equations are then solved for the fields as functions of T, /i„ and ^d- They give, 



dU_ 
dU_ 



(Pp nif 



f=u,d 



6ty: 



f=u,d 



(27r)3 
(2^ 



f = u,d 



Af{Ef)t 



-(Ef+f,s)/T 



Finally we note that the values of the parameters used are exactly the same as used in Ref. [30] 



(Al8a) 
(A18b) 

(A18c) 

(AlSd) 
(AlSe) 
(A18f) 
(A18g) 
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